Exact solution for core-collapsed isothermal star clusters 
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Star clusters in isothermal spheres are studied from a thermodynamic point of view. New density 
profiles are presented, that describe the collapsed phase at low temperatures. At the transition a 
set of binaries is formed that carries 13 % of the gravitational energy, while also a huge latent 
heat is generated. The total energy of the binaries is fixed by thermodynamics. In the canonical 
ensemble all specific heats are positive, while previously discussed negative specific heat solutions 
are metastable. In a microcanonical ensemble the latter remain partly dominant, 
keywords: ISM: globules, stellar dynamics, gravitation 
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Globular star clusters typically have 10 5 -10 6 stars of 
solar mass. In our Galaxy about 150 of them have been 
observed, while their total number is estimated to be 
several hundreds. Their most intriguing property is core 
collaps: Some 20 % of the globulars in our Galaxy have 
undergone a "violent relaxation" towards a collapsed 
state jjj . The core collaps is stopped by the formation of 
strongly bound binary stars ("hard binaries"), and possi- 
bly also more exotic objects like blue stragglers or black 
holes, that provide energy to the remaining core || 

Self-gravitating star clusters confined to isothermal 
spheres have served as important models to understand 
the evolution of globular clusters. The problem was stud- 
ied extensively by Chandrasekhar Q] . Todays status was 
presented in an overview by Padmanabhan ||. When 
varying the density contrast between center and bound- 
ary, do — n(0)/n(R), there occur three regions: i) the be- 
havior is regular for do < 32.2; it) for 32.2 < do < 702.9 
the solutions have a negative heat capacity; Hi) for 
do > 702.9 the solutions are unstable j|] Q. 

The occurrence of negative heat capacities is counter- 
intuitive, and has led to violent discussions. Let us stress 
that we are discussing an equilibrium situation, which, 
in principle, should behave normally. Out of equilib- 
rium the situation would be different. Negative heat 
capacities do occur for black holes ||] and even for sim- 
ple model glasses f|, without any contradiction of the 
non- equilibrium thermodynamics (§] ||. In a canoni- 
cal equilibrium description the heat capacity is propor- 
tional to the energy fluctuations, and necessarily posi- 
tive. Thirring pointed out that this objection does not 
apply in a microcanonical approach ffij. Lynden-Bell 
and Lynden-Bell considered the problem of negative heat 
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situation has nevertheless remained unsatisfactory. 

As a standard thermodynamic description could be 
given for black holes ||, we have started to believe that 
thermodynamics applies to gravitation and cosmology in 
general. Here we will set as our goal the description of the 
isothermal star cluster from a standard thermodynamic 
point of view. By this we mean: consider the lowest free 



energy, the specific heat, etc. when adiabatically lowering 
the temperature, thus employing a canonical ensemble. 
The number N ^> 1 of identical stars, the total mass 
M = Nm and the sphere radius R are kept fixed. 
The temperature and energy scales are set by 
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Let us introduce the reduced temperature 
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The lowest temperature for which solutions are known 
is T c = 0.396741 Tq g, corresponding to r c = 0.79346. 
We shall address the open question how the system 
behaves at lower T. It should be described by the 
Poisson equation for the gravitational potential cj>(r) = 
—Gm J d 3 r' n(r')/\r — r'|. For isothermal spheres one has 
n(r) = Ae-P m ^ r \ leading to 



V 2 </>(r) = 4irGmn(r) = AirGmAe 



-/Jm^(r) 
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where n is the number density and (3 = 1/T. The regular 
solutions are spherically symmetric and can be mapped 
on the Emden equation 
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with y(0) = y'(0) = 0. They are well understood, see ||. 

Chandrasekhar Q and later Lynden-Bell and Wood 
mention the singular solution 
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Pm<j){r) = 2 In 2. (5) 
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Though these authors do not pay much attention to it, 
it will become the corner stone of our analysis. The so- 
lution (H) has t = 1, and thus only exists at T = Tq/2. 
To describe the physics at low T, we have asked our- 
selves: Can eq. (||) be deformed continuously for t^I 
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? We have therefore considered singular solutions with 
only cylindrical symmetry 
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subject to the north-south symmetry 8 — > n — 6, imply- 
ing h(c) = h(—c). We call this the "core" solution, as 
there is no way to identify a halo in our static sphere 
with identical masses. The gravitational potential reads 



f3m<f>{r, I 



2 In— + h{cos$)--- \ dch(c) (7) 



where r is to be fixed by the form of the solution. With- 
out loss of generality we may require J^dch(c) = 0. 
From (||) it follows that h has to obey 



2 + (1 - c 2 )h"{c) - 2cti(c) = -e- 1 " e-' i( 
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with Iq = In f Q dce~ h . As it should, the case h = leads 
us back to eq. (Q) and r = 1. 

It is advantageous to introduce a new variable by set- 
ting c = tanhi, and a new function x(t), by defining 



h(c) = x(t) +ln(l - c 2 ) 
Some algebra shows that eq. (H) simplifies to 
d 2 x 2 



dt 2 
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This is just the equation for one-dimensional motion in 
an exponential potential. The general solution is 



x = 21ncosh[6(f - t )] - ln(6 2 r) - 7 



(11) 



Symmetry around c = imposes = 0. The parameter 
6 follows from the behavior for t — ► oo or c — * 1. The 
latter can be determined by integrating ( ra) from to c, 
using /i'(0) = 0. Comparing with @ and ( jll| ) then yields 
6 = 1/t. The solution of eq. (@) thus reads 
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The constant follows finally by imposing dch = 0, 
yielding 



7 = 2A4(-)-2 + lnr 

T 

The normalization of the /i-function (|l2|), 
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cannot be expressed in elementary functions. It has 
the behavior Nh(v) = In 2 + (7r 2 /24)^ 2 for small ^; 
/Vk(l) = 1, A/£(l) = |, while A/^) -» ^ln2 + 7r 2 /(24iy) 
for large ^. From the second equality it is easily seen 
that v 3 M^(v) = Af{l(l/v), which implies the symmetry 
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It is important to know when the solution is stable. 
Let us consider angular variations of the form n(r) — ► 
n(r, 6)[1 + f(y)], where y — atanh(cos 0)/t, and satisfy- 
ing J d 3 r nf = 0. For small / the free energy is invariant 
to first order, while to second order it reads 
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with n' = n(r' ,9') and /' = f(y'). For our profile the 
term between brackets leads to the Schrodinger equation 



f"(.v)--%z-f(v)=Ef(v) 
cosn y 
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For r > 1 there is a bound state, / = Asmhy/ cosh A y 
with A = (-1 + Vl + 8r)/2. Since £ < such fluc- 
tuations would lower the free energy, thus making the 
solution unstable. This has a simple physical interpre- 
tation: if r > 1 the density diverges on the north-south 
axis; this unnatural situation is now seen to be unstable. 
If no other instabilities occur, the solution (||), jl2]) is 
thus stable for < t < 1. 

Let us now see what is the physical meaning of the 
new solution. When coming from high T the system is, 
of course, in the regular solution. It can be described in 
terms of the Milne variables (u, v), with v = (3Tq — 2/t, 
that satisfy gj § 
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with initial condition u(0) = 3. In the u-v plane the so- 
lution spirals counter clockwise around the point (1,2), 
corresponding to infinite density contrast do- 
The entropy is generally defined as 
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where A(T) ~ T 1 ' 2 is the thermal wavelength. We shall 
subtract the constant S = N \u[V2t: R 3 / N A 3 (T G )] from 
S in order to condense our notation. 

The high temperature phase then has the properties 
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We can now identify the binodal temperature Tb - ,, 
and the spinodal temperature T$ — T c . In the interval 
(Ts,Tb) (i.e. r c < t < 1) there are two stable phases: 
the regular one and the new core state. To lower the free 
energy, the system may collaps at Tg. The free and in- 
ternal energy, the entropy and the specific heat are then 
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When we lower T the system can also stay in the under- 
cooled regular state, and undergo a core collaps at some 
transition point T t — \tiTq between T$ and Tb- If it 
did not yet do so, the system must collaps at Tg, simply 
because there is no other option. 

At each of these transition points the free energy is 
discontinuous, a truly uncomfortable feature. The solu- 
tion to this paradox is to assume that also binary stars 
are formed with free energy bridging this gap. The col- 
laps can thus be seen as a phase separation of the reg- 
ular phase into a core and binaries. It is well known 
that only a few hard binaries can be involved. Indeed, 
a finite central mass fraction would lead to a behavior 
[3m(j) ~ —R/r, which is incompatible with the Poisson 
equation (||). We therefore assume that the entropy of 
the binaries can be neglected, even though they must 
have a large energy. Using J7binaries=-Fbmaries we obtain 
from the difference between (poj) and (p4|) 
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where ut is the u- value corresponding to v = 2/t{. The 
above expression ranges from —0.0765 at Tb to —0.10193 
at Ts, and is roughly equal to 13 % of the pre-collaps 
gravitational energy. 

As usual at first order phase transitions, there is also 
a latent heat 

^Q = l + (ut -3 +Af h (T t )-l In utr^n (29) 
UG 2 



At Tb this equals 0.4062, being no less than 60.6 % of 
the precollaps gravitational energy, while at Ts it equals 
0.22043, still 27.8 %. These striking effects are induced 
by the strong binding nature of the core, expressed by its 
1/r 2 density profile (pf). It leads to a low value of the core 
energy (|25|). This is the manifestation of Lynden-Bcll's 
violent relaxation |lj], and it is violent indeed ! 
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FIG. 1. Free energy F as function of temperature T. 
Below Tb the system may enter the collapsed phase, which 
has lowest free energy. The regular phase may also be un- 
dercooled; then collaps occurs at a lower temperature, and 
the modified binary energy shifts the collaps curve upwards. 
For the lowest possible collaps temperature, Ts, the collaps 
free energy is presented by the upper curve. The upper 
branch of the regular curve has a negative specific heat, and 
is metastable. 

From eqs. ( p6| ) and ( p7| ) it is clear that C core ^ 
TdS covc /dT, which looks quite disturbing. However, we 
should not forget that the binary content is also a thermal 
part of the system. Assuming that equilibrium thermo- 
dynamics applies to the full system, we infer 
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Neglecting ^binaries and using eqs. (g6j), (|27j), this yields 
= ^(r) + (l-rV^(r) (31) 
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It is positive, implying that upon lowering T more and 
more energy is stored in binaries. Together with ( p8| ) this 
shows that thermodynamics relates the energy of the bi- 
naries to the properties of the core. Both contributions 
to the specific heat are positive, and so is their sum 
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At T B one has C coro = N/2, Canaries = N, while at 
T s : Ccore = 0.79504 7V, C binari c S = 0.99432 N. Finally, 
at T = we find the Dulong-Petit type behaviors C core 
= (5/2 - In 2)iV, C binarics = In 2 A. 

The total internal energy of the collapsed phase reads 
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The final effect of metastability, leading to a T t < Tg, 
is to make the energy content of the binaries smaller. In 
other words, knowledge of their energy allows the deter- 
mination of the collaps temperature T t . 

In Figure 1 we plot the total free energies of the prob- 
lem. Our approach solves the often discussed paradox 
of negative heat capacities of certain regular solutions. 
Indeed, as shown in Fig. 1, for T in a region above 
Tb there are two regular solutions. The one with larger 
free energy is also the one with the largest density con- 
trast (do > 32.2), and a negative heat capacity. It is 
metastable and cannot be reached by adiabatic cooling. 

In Figure 2 we present the specific heat as function of 
temperature for the regular solution, as well as for the 
binaries, the core and their sum. 
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FIG. 2. Specific heat C as function of temperature T. All 
contributions are positive. Upon cooling the system may col- 
laps at or below Tb- The largest chance occurs at Ts, where 
the regular solution has diverging energy fluctuations. 



capacity, but that is not forbidden in the microcanonical 
ensemble. For the energy of the collapsed phase we must 
adjust onset value for the binaries, leading to the binodal 
value Ub = —0.2679 Uq- The spinodal point is the point 
of instability of the regular solution, Us — —0.3342 Uq- 
The microcanonical ensemble thus exploits the regular 
solutions more. But, as seen from Figure 3, as soon as 
the collapsed solutions are stable, they are also in the 
microcanonical approach the dominant ones. 

In conclusion, we have presented new stable density 
profiles that describe the core collapsed phase of isother- 
mal clusters. Using standard thermodynamics the energy 
content stored in binaries is determined. In the canonical 
ensemble all specific heats are positive. Standard notions 
such as binodal and spinodal temperatures, phase sepa- 
ration, undercooling and latent heat find their natural 
place, just as in condensed matter physics. 

Our results probably pertain to realistic clusters, 
whose core is to a good approximation isothermal. Also 
there we expect that thermodynamics fixes the energy of 
the binaries. For clusters with north-south axis perpen- 
dicular to our line of sight, the projected density of the 
profile @ is proportional to e~' l ( cos e ) acos(r/ R) / V . Such 
a 1/r law is known from observations jl2] ] J3j . 
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FIG. 3. Entropy S as function of energy U . Coming from 
high U, the system enters the region with negative heat ca- 
pacity at U c = —0.1965 Ug- At Ub it can go to the collapsed 
state. If it fails to find it, it can ultimately collaps at Us, 
where the regular solution becomes unstable. 
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Let us finally discuss the microcanonical situation. We 
have to consider the entropy curves S(U) on lowering U. 
For the regular solution it follows by eliminating T be- 
tween ( p2| ) and (^lj) . The entropy of the stable collapsed 
phase has a maximum S^uaps^ = 1) = \- The reg- 
ular phase has at that entropy already a negative heat 
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